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Optical spectroscopy study of three dimensional Dirac semimetal ZrTe 5 
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Three dimensional (3D) topological Dirac materials are under intensive study recently. The layered compound 
ZrTe 5 has been suggested to be one of them by transport and ARPES experiments. Here, we perform infrared 
reflectivity measurement to investigate the underlying physics of this material. The derived optical conductivity 
exhibits linear increasing with frequency below normal interband transitions, which provides the first optical 
spectroscopic proof of a 3D Dirac semimetal. Apart from that, the plasma edge shifts dramatically to lower 
energy upon temperature cooling, which might be associated with the consequence shrinking of the lattice 
parameters. In addition, an extremely sharp peak shows up in the frequency dependent optical conductivity, 
indicating the presence of a Van Hove singularity in the joint density of state. 

PACS numbers: 71.55.Ak, 78.20.-e,72.15.Eb 


Three dimensional (3D) topological Dirac materials, such 
as Dirac semimetals and Weyl semimetals, have attracted 
tremendous attention in condensed matter physics and mate¬ 
rials science in recent years. The 3D Dirac semimetals pos¬ 
sess 3D Dirac nodes where the valence and conduction bands 
touch along certain symmetric axis. The 3D Dirac node is 
protected against gap formation by crystalline symmetry and 
can be considered as two overlapped Weyl nodes of opposite 
chirality [1, 2]. The materials are expected to host many in¬ 
triguing novel phenomena and properties, such as the appear¬ 
ance of Fermi arcs on the surfaces [3], giant diamagnetism 
[4, 5], and very large linear magnetoresistance [6, 7]. When 
time reversal symmetry is broken, a 3D Dirac semimetal can 
be driven to a Weyl semimetal. This has stimulated strong 
interest to explore chiral magnetic effect [8-15], i. e. the 
charge-pumping between Weyl nodes of opposite chiralities 
under magnetic field, which is a macroscopic manifestation 
of the quantum anomaly in relativistic field theory of chiral 
fermions. 

Only a few materials were predicted to be 3D Dirac 
semimetals. Among them, NajBi and Cd 3 As 2 are the two 
best known compounds [16, 17]. They were confirmed to be 
3D Dirac semimetals by angle-resolved photoemission spec¬ 
troscopy (ARPES) experiments [18-21]. Very recently, an¬ 
other layered compound, ZrTes, was suggested to be also a 
3D Dirac semimetal based on transport and ARPES experi¬ 
ments [13]. ZrTe 5 crystallizes in the layered orthorhombic 
crystal structure, with prismatic ZrTe () chains running along 
the crystallographic a-axis and linked along the c-axis via 
zigzag chains of Te atoms to form two-dimensional (2D) lay¬ 
ers. Those layers stack along the b-axis. Most attractively, the 
chiral magnetic effect was clearly observed for the first time 
on ZrTes through a magneto-transport measurement [13]. 

The 3D band dispersions in Dirac semimetals enable the 
bulk measurement techniques, e.g. the optical spectroscopy 
measurement, to provide key information about the charge 
properties of the materials. For any electronic material, such 
information is of crucial importance for understanding under¬ 


lying physics. The situation would be different from the in¬ 
vestigations on 3D topological insulators where the focus is on 
the non-trivial surface state. For those materials, the bulk mea¬ 
surement technique could hardly detect or separate the con¬ 
tribution from the 2D Dirac fermions of topological surface 
states, as a result, the experimental probes were largely lim¬ 
ited to surface sensitive measurements like ARPES or scan¬ 
ning tunneling microscopy (STM). 

Although a 3D Dirac semimetal could be viewed as a bulk 
analogue of 2D graphene in terms of band dispersions, a 
key difference exists in charge dynamical properties between 
2D and 3D Dirac fermions. For 2D Dirac fermions like in 
graphene, the real part of conductivity is well known to be 
a frequency independent quantity [22], while for 3D Dirac 
fermions with linear band dispersions in all three momentum 
directions, the real part of conductivity grows linearly with 
the frequency [23, 24]. To our knowledge, no optical spectro¬ 
scopic measurement has been reported on 3D Dirac semimet¬ 
als. In this work we present temperature and frequency depen¬ 
dent optical spectroscopy study on ZrTes. Our study reveals 
several peculiar features. The density of free carriers keeps 
dropping with temperature decreasing and an unusual sharp 



FIG. 1: Temperature dependent resistivity of ZrTej along a-axis. 
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FIG. 2: Temperature dependent optical reflectivity R(oj) below 600 cm 1 along a- (panel a) and c-axis (panel b). Panel (c) displays R(u>) below 
10000 cnr 1 along a axis, while the inset shows the spectrum up to 50000 cnr 1 at room temperature. 


peak shows up in the reflectivity and conductivity spectra in 
the mid-infrared region. Of most importance, the linear ris¬ 
ing of the optical conductivity is observed in a relatively large 
energy scale, just as expected for a 3D Dirac semimetal. 

Single crystals ZrTe s were grown by using Te flux 
method. 100 grams high purity 7N (99.99999%) Te and 6 N 
(99.9999%) Zr were loaded into 20 mm diameter double- 
walled quartz ampoules with 200 mm length and sealed under 
vacuum. The composition of the Zr-Te solution for the ZrTes 
crystal growth is Zr 0 . 0025 Te 0 . 9975 , and the largest ZrTes single 
crystals obtained are around ~ 1 x 1 x 20 mm 3 . Detailed 
growth procedure was described elsewhere [13]. 

Figure 1 displays the dc resistivity in the chain direction 
as a function of temperature of ZrTes. The resistivity curve 
shows rather unusual features. With decreasing temperature 
from 300 K, p(T) decreases slightly first, then increases be¬ 
low roughly 200 K and reaches a peak near 60 K. The over¬ 
all small values of p(T ) signal the (semi)metallic nature of 
the compound, which are in consistent with the optical mea¬ 
surements presented below. In earlier reports, the resistivity 
peak located near 150 - 170 K [25, 26]. The much lower peak 
temperature could be attributed to a much lower defects in 
the present sample. The unusual features are likely associated 
with its characteristic electronic structure, which could be sen¬ 
sitive to the temperature-induced structural change. Earlier 
band structure calculations suggested presence of extremely 
small and light ellipsoidal Fermi surfaces, centered at the cen¬ 
ter (T point) of the bulk Brillouin zone (BBZ) [27], and quan¬ 
tum oscillation measurements indicated three tiny but finite 
Fermi surfaces [26], with the effective mass in the chain di¬ 
rection (m*a-0.03m ( ,) being comparable to that in a prototyp¬ 
ical 3D Dirac semimetal, Cd 3 As 2 [28]. More recent ab initio 
calculations indicated that ZrTes compound locates close to 
the phase boundary between weak and strong topological in¬ 
sulators [29]. Nevertheless, very recent transport and ARPES 
experiments identify it to be a 3D Dirac semimetal [13]. 

The polarized optical reflectance measurements with E//a- 
axis and E//oaxis were performed respectively on a combi¬ 
nation of Brake 113V, 80V and grating-type spectrometer in 
the frequency range from 30-50 000 cm -1 . The reflectivity 
along a- and c-axis in the far-infrared region are displayed in 


Fig. 2 (a) and (b), respectively. The a-axis reflectivity over 
broad energy scales is shown in Fig. 2 (c). There are several 
prominent phonon absorption features lying below 200 cm -1 , 
which get more pronounced upon cooling. At room temper¬ 
ature, the reflectivity R(a>) shows a well-defined plasma edge 
below 350 cm " 1 and approaches unity at low frequency. This 
indicates clearly a metallic response, which is consistent with 
the results of transport measurements. The edge frequency, 
being usually referred to as ’’screened” plasma frequency, is 
related to the density n and effective mass m* of free carri¬ 
ers by u)'p oc n/m*. The unusual low value of this frequency 
suggests a very small n. As the temperature decreases, the 
edge shifts to lower frequency, indicating a further reduction 
of the number of free carriers. Below 80 K, the extremely 
strong phonon signals tend to blur the plasma edge. In the 
meantime, the carrier scattering rate is also reduced as the re¬ 
flectance edge gets steeper. The combination of the two ef¬ 
fects is likely the driving force of the broad peak appeared at 
60 K in p(T). Furthermore, the rapid reduction of the plasma 
frequency yields optical evidence that the compound tend to 
approach the semimetal state at very low temperature. 

R(w) along the c axis shows similar behaviors as along the 
a axis. The plasma edge almost lies at the same position as 
along the a axis at different temperatures. However, the over¬ 
all reflectivity is a little bit lower along the c axis, which is 
absolutely reasonable because the conductivity along this di¬ 
rection is supposed to be lower. Moreover, the phonon signals 
shown in Figure 2 (a) and (b) are quite distinct, revealing dif¬ 
ferent oscillating modes. In the following paragraphes, we 
will only discuss the characteristic features along the a axis. 

Basically, the spectra at higher energies are dominated by 
interband transitions, as shown in Fig. 2 (c). The peaks cen¬ 
tered around 2700 cm" 1 ,4500 cm " 1 and 7800 cm " 1 represent 
for different transitions. Remarkably, the interband transitions 
exhibit temperature dependence. This is seen more clearly 
for the second peak. As temperature decreases, the peak fre¬ 
quency moves to higher energy and its spectral weight is en¬ 
hanced simultaneously. The enhancement is linked to the re¬ 
duction of Drude spectral weight. As we shall explain below, 
it may be caused by the shift of chemical potential associated 
with temperature-induced subtle crystal structural change. 
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FIG. 3: The optical conductivity at different temperatures up to (a) 
300 cm" 1 and (b) 10 000 cnr 1 . The short dashed lines are low fre¬ 
quency extrapolations by the Hagen-Rubens relation. The inset of 
the upper panel shows the temperature dependent evolution of the 
plasma frequency oj p and scattering rate of free carriers y by black 
solid squares and blue open squares respectively. 


The real part of optical conductivities cr \ ( oj) are derived 
from R(cn) though Kramers-Kronig relation, as displayed in 
Figure 3. The Hagen-Rubens relation was used for the low 
frequency extrapolation, and the x-ray atomic scattering func¬ 
tions were used in the high frequency extrapolation [30]. Fig¬ 
ure 3 (a) shows cri(cj) in the low frequency. Although the 
screened plasma edge can be clearly observed in R(cn), the 
Drude component in cri (oj) appears mainly in the extrapolated 
region. Apparently, the Drude spectral weight reduces signif¬ 
icantly at low temperature, being consistent with the dramatic 
shift of plasma edge towards low frequency in R(co). In ad¬ 
dition, three sharp phonon peaks are resolved at 46, 86 and 
185 cm -1 , respectively. When the temperature decreases, the 
phonon frequencies shift slightly towards higher energies. 

The free carrier contributions could be quantitatively es¬ 
timated from the application of Drude model, which has a 
peak centered at zero frequency and the peak width is just the 
carrier scattering rate y. The Drude spectral weight, which 
gives the plasma frequency oj p = \f4jmJlrF, could be alter¬ 
natively obtained from integrating the low frequency spectral 
weight by or p — 8 ^ c cr\doj, where u> c is the cut off fre¬ 
quency. The integration should sum up all the free carrier 


contributions, without the influence of any interband transi¬ 
tions. We chose the frequency where <xj (oj) reaches a mini¬ 
mum as a> c , which was identified to be 30, 30, 70, 130, 220 
and 330 cm" 1 in sequence corresponding to increasing tem¬ 
peratures. The yielded temperature dependent scattering rate 
and plasma frequency are plotted in the inset of 3 (a). The 
scattering rate decreases from 40 cm" 1 at 300 K to 4 cm" 1 at 
8 K. Meanwhile, a> p drops monotonically upon cooling from 
1400 cm" 1 to 300 cm" 1 . Assuming the effective mass of free 
carriers remains unchanged, then 95% of the population is re¬ 
duced. The dramatic reduction of free carriers could hardly 
been explained by reduced thermal excitations. Here, we pro¬ 
pose that this peculiar phenomenon is related to the chemical 
potential, which moves closer to the Dirac node as tempera¬ 
ture decreasing, along with the lattice shrinking. Theoretical 
calculation suggested that the electronic properties of ZrTes 
are extraordinarily sensitive to the variation of the lattice pa¬ 
rameters [29], 

Figure 3 (b) presents the optical conductivity in the energy 
scale up to 10 000 cm" 1 . Corresponding to R(«), an excep¬ 
tional sharp peak shows up at around 4000 cm" 1 , along with 
a shoulder-like feature at around 2800 cm" 1 . When tempera¬ 
ture decreases, the central frequency of the sharp peak shifts 
to higher energy and the half with of it gets smaller. At the 
lowest temperature, it even becomes quasi-divergence. This 
feature has never been observed in any other 3D materials. 
We noticed that the reflectance of single-walled carbon nan¬ 
otubes, whose joint density of states (JDOS) are constituted 
by individual Van Hove singularities, was reported to show 
geometry-dependent resonant peaks [31], which are identified 
to stem from interband transitions between these singularities. 
Consequently, it is most likely that the quasi-divergent peak 
exhibited in cri(to) also requires for a Van Hove singularity in 
the IDOS of ZrTes, which is rarely seen in a 3D bulk material. 
Nevertheless, the ZrTes single crystals are easily cleaved both 
along the b and c axis, forming a needle-like compound. The 
lattice parameters of a=3.9876 A, b=14.505 A, c=13.2727 A 
also indicate a quasi-one dimensional property, which tends to 
strengthen the effect of a Van Hove singularity. 

Usually Lorentz model is employed to describe the inter¬ 
band transition. However, for present compound, the conduc¬ 
tivity data above the Drude component could not be repro¬ 
duced by the Lorentz model due to a linear increase of con¬ 
ductivity with frequency. To make it clearly to see, the optical 
conductivity at 8 K is displayed in Fig. 4, where the linear 
increasing is fitted by a red dotted line. Baasi and Virosztek 
showed that, for a noninteracting electron system consisting of 
two symmetric energy bands touching each other at the Fermi 
level whose Hamiltonian depends on the magnitude of the mo¬ 
mentum through a power-law behavior, the optical conductiv¬ 
ity of the interband contribution has a power-law frequency 
dependence with cri(oS) oc )~ , where d is the dimen¬ 
sion of the system and z represents for the power-law term 
of the band dispersion [32]. For example, the two dimen¬ 
sional graphene with linear dispersion of e(k) oc |k| (that is, 
d-2 and z= 1) gives rise to a frequency independent conduc- 
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FIG. 4: The optical conductivity at 8 K below 1200 cm 1 . The red 
dotted line are the linear fitting of cr l (a>) 


tivity, which has been confirmed by optical experiments [22]. 
For the case of 3D compound ZrTes, the linear rising or (oS) 
demands for — - 1. By applying d- 3, a linear dispersion 
z- 1 is obtained, in good agreement with a Dirac semimetallic 
behavior. 

As a matter of fact, the linearly growing optical conduc¬ 
tivity has also been observed in a Rh-doped NdilriO? below 
50 meV [33]. Since pyrochlore iridates were suggested to 
host Weyl semimetals [34], the observation has been taken 
as the evidence of a Weyl semimetal even though the lin¬ 
ear increasing only persistd in a very small energy scale. 
Significantly, several quasicrystals (such as Al 63 .sCu 24 .sFei 2 , 
Al 75 . 5 Mn 2 o. 5 Siio.i) has shown similar behaviors as well [35], 
all of which are lack of inversion symmetries and were sug¬ 
gested to be candidates of Weyl semimetals. Furthermore, 
3D massless electrons with huge Fermi velocity are certified 
by linear rising optical conductivity in the zinc-blende crys¬ 
tal Hgi_ x Cd v Te, with a proper doping concentration close to 
x—0 .17 [36]. What we observed here, however, is the first bulk 
evidence of a time reversal and space inversion symmetry pro¬ 
tected 3D Dirac semimetal in the as-grown single crystal. 

The linear increase of optical conductivity is an unparal¬ 
leled property of bulk 3D Dirac band dispersions. Theo¬ 
retical studies show that the optical conductivity contributed 
by per Weyl node (lifting the spin degeneracy) is [23, 24] 
ci ( w ) = f^f7 = ^°2TT’ where = 2e 2 //t=7.727xl0~ 5 Q~ 1 is 
the quantum conductance. For ZrTes, band structural calcula¬ 
tions indicates only one possible Dirac node near I point per 
ZrTes layer [29]. Considering that each Dirac node has a spin 
degeneracy of 2, we can get the Fermi velocity vf — 2.59xl0 6 
cm/s in terms of the slope of the conductivity shown in Fig. 
4. This value is significantly smaller than the Fermi velocity 
determined by ARPES experiment along a- or c-axis (smaller 
but close to 10 8 cm/s). We remark that, for 3D Dirac semimet¬ 
als, although the linear dispersions exist for all three momen¬ 
tum directions, their slopes or Fermi velocities are in fact dif¬ 
ferent. ZrTes has a layered structure, the Fermi velocity in 
the direction perpendicular to the layers (along b- axis) must 


be much smaller than the other two directions, unfortunately 
presently available measurement has not yielded such infor¬ 
mation. The Fermi velocity determined from optical conduc¬ 
tivity for a 3D Dirac semimetal must contain contributions 
from all momentum directions. 

For a rigorous Dirac semiemtal, the Fermi energy lies right 
across the Dirac points and the Fermi surfaces would shrink 
to isolated points in the BZ. As a result, there would be no 
free carriers. The optical conductivity of ZrTes shows obvi¬ 
ous Drude component at high temperatures, implying that the 
Fermi energy stays away from the Dirac point. With temper¬ 
ature decreasing, the density of free carriers decreases mono- 
tonically and almost all of them are lost at the lowest temper¬ 
ature, which is consistent with the approaching of the Fermi 
level to the Dirac point. ARPES measurement has produced a 
very similar scenario. It is also worth noting that the frequency 
where the conductivity deviates from linear dependence can 
not be treated as the cutofF energy of linear band dispersion, 
as ct\(u)) is dramatically modified by other interband transi¬ 
tions at higher energies. 

In conclusion, we have studied the optical conductivity of 
the 3D topological Dirac semimetal material ZeTes. The num¬ 
ber of free carriers are proved to be much smaller than con¬ 
ventional metals, manifesting a semimetallic behavior. Fur¬ 
thermore, it decreases as temperature decreasing, probably 
caused by the modification of chemical potential, which might 
be sensitive to the lattice parameter variation. An unexpected 
quasi-divergentpeak is also observed in the optical conductiv¬ 
ity, which is possibly associated with a Van Hove singularity 
in the JDOS. Most importantly, the optical conductivity grows 
linearly with frequency below 1200 cm -1 , presenting the first 
optical spectroscopic evidence of a 3D Dirac semimetal. 
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